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Abstract 



Experiments on surface melting of several organic materials have shown 
contradictory results. We study the Van der Waals interactions between in- 
terfaces in surface melting of the bulk CH4 and interfacial melting of the CH4 
film on the MgO substrate. This analysis is based on the theory of Dzyaloshin- 
skii, Lifshitz, and Pitaevskii for dispersion forces in materials characterized by 
the frequency dependent dielectric functions. These functions for magnesium 
oxide and methane are obtained from optical data using an oscillator model of 
the dielectric response. The results show that a repulsive interaction between 
the solid-liquid and liquid-vapor interfaces exists for the bulk methane. We 
also found that the van der Waals forces between two solid-liquid interfaces 
are attractive for the film on the MgO substrate. This implies that the 
van der Waals forces induce the presence of complete surface melting for the 
bulk methane and the absence of interfacial melting for CH4 on the MgO 
substrate. 
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I. INTRODUCTION 

As the temperature is increased from below the bulk triple point, some solids begin to 
melt from their surfaces inward. If the thickness of the melted layer increases without limit 
as the temperature approaches the bulk triple point, then the solid is said to undergo surface 
melting The melting is incomplete if the thickness of the melted film remains finite as 
the temperature increases at all temperatures including the triple point. A special case of 
incomplete surface melting where the liquid film thickness is zero is sometimes referred to as 
surface nonmelting. For example, it was found that the P6(110) surface melts. The P6(lll) 



surface does not melt, whereas P6(100) exhibits incomplete melting. Therefore, incomplete 
or blocked melting appears on (100) fee metal surfaces with an intermediate packing density 
between (110) and (111) characterized by a quasi-liquid film of the finite and 

small thickness. Surface melting was studied for several metals and other materials, such 
as ice, germanium, and multilayers of gases adsorbed on graphite |]1],0. High-temperature 
disordering for surfaces of several organic materials was also studied. 

X-ray reflectivity measurements were used for high-temperature studies of the (100) face 
of caprolactam {CqHuON) 0,0. Initially a series of layering transition followed by prewet- 
ting was suggested, but later a surface hydration mechanism was proposed for temperature 
disordering of this surface. Chernov and Yakovlev studied surface melting for the (001) and 
(010) faces of biphenyl (CuHiq) using the ellipsometry technique 0. They found that the 
thickness of a quasi-liquid layer at 0.5 K below the bulk melting point reaches lOOA for 
the (001) face and 300A for the (010) face. In contrast, de Jeu and coworkers using the 
synchrotron X-ray reflectivity measurements have shown the absence of surface melting for 
the (001) face of the bulk biphenyl up to 2 K below the melting point |Q. It was proposed 
that surface melting observed by Chernov and Yakovlev is induced by the roughness of the 
glass used in their ellipsometry experiment ||^. Surface properties of biphenyl, caprolactam, 
and other organic materials and organic thin films at high temperatures deserve further 
experimental and theoretical studies. Investigations of organic films are also important for 
applications, such as tribology and lubrication, corrosion inhibition, modification of the sub- 
strate optical properties, molecular electronics, and various industrial chemical processes. 

Surface melting for a ~ 10 layer thick film of CH^ adsorbed on MgO{100) was studied by 
neutron scattering P,p!0| . It was found that the thickness of the quasi-liquid phase changes 



from one layer at 72 K to ~ 6 layers at 90.3 K. (At atmospheric pressure the bulk melting 
point of methane is 90.66 K, whereas the liquid-vapour transition occurs at 111.56 K.) The 
neutron-scattering measurements of the translational diffusion coefficient above 72 K have 
shown that this quantity is in the 10~^ cm^ range, i.e., it has a value typical for mobility 
in quasi-liquid films. It was found that for a methane thick film (~ 20 layers) at T = 94 K 
(i.e., above the bulk melting temperature) almost all layers are liquefied 0]. Recently the 
layering and melting properties of methane on MgO{100) were studied between 70 K and 96 
K using high-resolution adsorption isotherms and for the film thickness of up to five layers. 



ri| . Melting transition was observed at 80 K for the monolayer, and at 85 K for layers 2 
through 4. 

Related to the experiments on surface melting of bulk caprolactam and biphenyl, it would 
be interesting to study surface melting of a bulk methane. These studies should determine 
whether surface melting of the CH4 film on MgO is inherent feature of this organic material, 
or it is induced by the substrate. Interfacial melting is an appearance of a liquid film at 
the interface between one solid material and the substrate of the other material as the 



bulk melting point of the first material is approached |l2| , p!3| . Therefore, surface melting is a 
special case of interfacial melting where the substrate is the vapor phase of the first material. 
In this work long-range interactions important for surface melting of a bulk methane, as well 
as for interfacial melting of a methane film on the MgO substrate are analyzed using the 
theory of Dzyaloshinskii, Lifshitz, and Pitaevskii (DLP) [Q. Our findings revealed that 
complete surface melting exists for the bulk methane. Because of an attractive interaction 
between interfaces melting is absent at the boundary between the CH4 film and the MgO 



2 



substrate. 

In the following the DLP theory of the van der Waals interaction between planar inter- 
faces and its consequences for surface melting are described in Sec. II. The calculation of 
frequency dependent dielectric functions of Cif4 and MgO is discussed in Sec. III. Results 
and discussion for the surface and interfacial melting problems are presented in Sec. IV. A 
summary and conclusions are given in Sec. V. 



II. SURFACE MELTING AND VAN DER WAALS FORCES 

The surface free energy per unit area for a system with a quasi-liquid film of the thickness 
/, at the temperature T is ||l|,3, P^] 

T) = + + L„(l - T/T^)l + + Fi,. (1) 

In this equation 7^/ and 7i„ are the free energies of the solid-liquid and liquid- vapor interfaces, 
respectively, is the latent heat of melting per unit volume, and is the bulk melting 
temperature. The fourth term Fg^ represents the gain in interfacial energy associated with 
the wetting of the solid surface by a quasi-liquid film. This term is the result of short- 
range interactions. Their contribution is oscillatory and decreases exponentially to zero 
for large /. Short-range forces produce a logarithmic temperature dependence of the liquid 
film thickness as the bulk melting temperature is approached. Long range van der Waals 
interactions are given by the last term Fi^, in Eq. (1). The contribution of the van der Waals 
forces for non-retarded interactions (where / is much smaller than the wavelength of light for 
a characteristic excitations) is given by Hl~^. The parameter H is the Hamaker constant 
and gives information about the magnitude of the van der Waals interactions in different 



media [|T^]. (Sometimes H' = 12-kH or H" = —12ttH are called the Hamaker constant.) 
Surface melting occurs when H > 0, i.e., when two interfaces (solid-liquid and liquid- vapour) 
repel each other. Then Fir{l) is monotonously decreasing function, or its minimum is not 
deep enough to bind the interface. The van der Waals contribution produces a power-law 
temperature dependence of the liquid film thickness as T ^ 

LmiTm — T) 



2HT, 



m 



For H < two interfaces attract each other and the melting is at most incomplete. Then a 
deep minimum exists on the free energy curve Fir, and the interface is trapped there. 

Dzyaloshinskii, Lifshitz, and Pitaevskii developed the theory of long-range dispersion 
interaction between macroscopic bodies ||l^JT6[] . The van der Waals interaction between 
interfaces is a consequence of electromagnetic fluctuations in a polarizable medium and the 
discontinuity of the dielectric function across the boundaries. In the DLP theory the force 
between interfaces (such as those shown in Fig. 1) is a function of the frequency-dependent 
dielectric functions e{uj) of materials in contact. In the following all media are taken as 
isotropic and continuous. The free energy per unit area of two media (1 and 2) separated 



by a film of thickness / (media 3, see Fig. 1) is |jT^,|T^,|T2 
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where 



and 



+/n(l - |^^"-^^"^|i^^"-^^"\^ exp(-x))), (3) 
(e2X + e3X2)(eiX + e3Xi) ^ ^ 



2lU{e,f'^/c, (4) 



^i,2=(x^-r^(l-^))'''. (5) 

The dielectric functions ei{uj), e2{uj), and e^ioj) are evaluated at imaginary frequencies i^m = 
i{27rkT/h)m. The functions are positive, real and decrease monotonically from the 
static dielectric function eo for ^ = 0, to 1 as ^ ^ oo. In the above expressions k, h, and 
c are the Boltzmann and Planck constants and velocity of light, respectively. The prime 
on the sum in Eq. (3) means that the term m = has to be multiplied with i. In the 
derivation of Eq. (3) retarded interactions were assumed. If retardation is ignored, then 
c ^ oo, — > 0, and Eq. (1) simplifies to 

rni T) = ^y Tdxxlnd- (^^(^^-) - ^Mm)){ed^U) - '^(^^rn)) \ 
^ir l^, ) 8^^2 l^J^ (,^(,^^) + e,{^U)){el{^U) + e^U)) V ' 

A common approach is to obtain an approximate value of e{uj) from the optical measurements 
T^ , pJ| -pT|. An oscillator model of the dielectric response is then used to generate data, i.e., 



the dielectric function is represented by 

where /j, Cj, and Qi are fitting parameters. The DLP theory was recently applied to surface 
melting of ice |T^,|T^, germanium, and several metals [^,0. Unfortunately, optical and 



dielectric properties of many organic materials are less studied than, for example of metals, 
semiconductors, or H2O. As a consequence, we were not able to find necessary optical 
data for biphenyl and caprolactam, two organic materials for which experiments on surface 
melting were performed. In contrast, we found data for optical properties of organic material 
CHi and dielectric substrate MgO p2|p3[ . 



III. DIELECTRIC FUNCTIONS OF METHANE AND MGO 



The complex dielectric functions e(z^) were constructed using experimental data for 
optical constants of CH^ and MgO [p2| , ^ . As suggested in Ref. p4| we minimized the 
objective function 



N 



exp / 
er [UJi 



+ 



1 2 



(8) 
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where e,.(u;) and ei{uj) are the real and imaginary parts of a dielectric function given by 
Eq. (7), and ^exp' denotes experimental values taken from Refs. |^,23|. Therefore, the 



calculation was performed by finding a minimum of an objective function based on a si- 
multaneous fitting of both real and imaginary parts of the dielectric function. We found 
that, in comparison with some other functions, the functional form in Eq. (8) is suitable 
for minimization. In principle, the most important factor in the minimization process is the 
quality of experimental data, i.e., the number of points and their distribution. It is known 
that the main contribution belongs to data in the UV region of the spectra pO|,^. The UV 



part of the spectrum for MgO is represented with many points P3. In contrast, in the UV 



spectra for the solid and liquid methane few points exist |22|. A smaller number of points 



gives more freedom in the choice of the parameters which minimize Eq. (8). Therefore, 
more local minima exist and they produce similar values of the objective function. This 
is a limitation for the calculation of dielectric properties of methane. Fortunately previous 
calculations for other materials have shown that the complex dielectric functions e{iC,) which 
enter expressions in the DLP theory are not very sensitive to the fitting procedure and the 
detailed experimental information. As suggested in Ref. ||2^ and followed by others |T7| , p!2 



we used the constraint er = n^isi where rims is the refractive index in the visible part of 
the spectrum. The value n^is is known with a great accuracy and fitting with this con- 
straint gives better results. Objective functions were minimized using the simplex method 
Il25| , |26[| . The accuracy of calculation was checked by changing minimization routines and 
also by calculating dielectric functions for other materials where extensive data are available 

Martonchik and Orton, motivated by astrophysical applications, recently calculated the 



optical constants of liquid and solid methane using the spectra from the literature [22 



Liquid methane was studied at the boiling and melting temperatures, whereas the solid 
phase I of methane was studied at the melting point and at T = 30 K. (The solid phase I is 
the face centered cubic phase of methane which at atmospheric pressure exists for T > 20 
K.) Therefore, data obtained by Martonchik and Orton for the solid and liquid methane at 
the melting point were used to model the solid and liquid phase. For the vapor phase we 
put e = 1. The static dielectric function of the gas methane is e = 1.000944 [^. The fits for 
CH^ were performed to match dielectric functions at zero frequency. These static dielectric 



functions are eo = 1.67 for the liquid methane, and eo = 1.74 for the solid methane p2 
The calculated complex dielectric functions of the solid and liquid methane are shown in 
Fig. 2. Corresponding fitting parameters are given in Tables I and II. The curves for 
the solid and liquid methane are close to each other, but do not cross. In similar calculation 
for surface melting of ice it was found that the complex dielectric functions of water and ice 
cross at 10 eV (i.e., for m ~ 205) 

The detailed review of optical properties of MgO was given by Roessler and Huffman 
p3| . The values of optical constants were tabulated in the range (0.002 — 586) eV. Because 
of a large set of data we use only selected points. From these data the function e{i^) for 
MgO was constructed. The static dielectric function of MgO is eo = 9.8 Here, as in 
Ref. where the fitting based on a double band was performed, we did not match the 
static dielectric function. The value of static dielectric function of MgO is much larger than 
the average value of in the UV part of the spectrum. Additional problem is the existence 
of a close sharp peak at Tiui = 7.8 eV. Therefore, the inclusion of the static dielectric function 
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into the fit for MgO prevents a good beliavior of the whole fitted curve. Fitting parameters 
for MgO are given in Table III. The complex dielectric functions e(«0 of MgO is shown in 
Fig. 2. Our four-band result is almost the same as a double band calculation from Ref. p8 . 



IV. SURFACE AND INTERFACIAL MELTING 



A. Bulk methane 



In surface melting problem the frequency dependent dielectric functions of solid, liquid 
(melt), and gas methane are respectively: ei = 6^(0;), es = ei{uj), and €2 = 1 (see Fig. 
1). The integral in Eq. (3) was calculated numerically using the Gaussian method. The 
accuracy of the integration technique was checked by changing numerical procedures and 
also by reproducing the results for ice |T^. The calculated free energy as a function of liquid 
film thickness is shown in Fig. 3. The minimum is reached for / — > 00. This functional 
form of F{1) is typical for complete surface melting. We found that the Hamaker constant 
is if = 5.26 X 10"^^ J. Similar values were calculated for metals with melted surfaces, such 
as Pb, AL and Au WEA. 



B. Methane film on MgO 



In interfacial melting problem, and for interfaces shown in Fig. 1: ei = eMgo{(^) is 
dielectric function for MgO, €2 = es{uj) is for a solid phase, and €3 = ei{uj) is for a liquid 
phase of methane. Therefore, we consider the solid Cif4 - liquid CH4 - MgO substrate 
system. The same numerical procedure as for the surface melting problem was applied. The 
free energy is shown in Fig. 4. The Hamaker constant is H = —6.03 x 10~^^ J. Therefore, 
the van der Waals forces produce an attractive interaction between interfaces and prevent 
interfacial melting. The functional form of F{1) shown in Fig. 4 was not found for interfacial 
melting of ice on various substrates [|T^], but two other forms of F{1) were calculated there. 
They correspond to complete and incomplete interfacial melting of ice on the particular 
substrate. 



V. SUMMARY AND CONCLUSIONS 



The interaction of interfaces present in surface melting of the bulk methane and interfa- 
cial melting of the CH4 film on the MgO substrate is analyzed within the Dzyaloshinskii, 
Lifshitz, and Pitaevskii theory of the van der Waals interaction. Frequency dependent di- 
electric functions for these materials are calculated from optical data. It was found that 
complete surface melting exists for methane. Our study is the first investigation of surface 
melting properties of the bulk methane. We also found that the van der Waals interactions 
do not produce interfacial melting of the CH^ film on the MgO substrate. Therefore, pre- 
melting observed in the experimental investigations of the CH4 films on the MgO substrate 
P^TTI is inherent feature of the bulk methane. Methane is till now the only organic material 
for which clear complete surface melting was found. The results obtained here depend on 
the quality of the optical measurements used to model dielectric functions. Nevertheless, 
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these results shed hght on the interaction between interfaces in surface and interfacial melt- 
ing involving methane. These results suggest experimental studies of surface melting of the 
bulk methane, for example using the synchrotron X-ray reflectivity. In addition, ab initio 
molecular dynamics study of the high-temperature properties of methane surfaces is feasible, 
similar as a recent simulation of a high-pressure behavior for bulk CH4 [^]. It is impor- 
tant to do all necessary optical measurements for biphenyl and caprolactam, and then the 
calculation of the van der Waals interactions for these materials. There is a possibility that 
surface melting observed for biphenyl on the glass surface is interfacial melting induced 
by the dielectric properties of the substrate. Till now the roughness of the glass and impurity 
induced interfacial melting were proposed to explain the results of Chernov and Yakovlev. 
A possibility of interfacial melting induced by the van der Waals forces for biphenyl on the 
glass also explains the absence of surface melting for the bulk biphenyl found using the X-ray 
reflectivity p|. Recent atomic force microscopy study of the (001) surfaces of n — C23H4S 
paraffin crystals \W\ shows that this technique gives a good opportunity to investigate melt 
growth on a molecular scale for surfaces of methane and other organic materials. 
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FIGURES 



FIG. 1. Schematic representation for configuration of interfaces between media characterized 
by the frequency dependent dielectric functions. 

FIG. 2. Complex dielectric functions. Dashed line for MgO is the fit from Rcf. [28]. The 
index m is defined in Eqs. (3) and (6) in the text, and m = 1000 corresponds to 48.8 eV. The 
contribution of various m-dependent terms in Eqs. (3) and (6) is obvious. 

FIG. 3. The free energy per unit area (oq is the Bohr radius) as a function of the film thickness I 
for surface melting of the bulk methane. The dashed line shows the free energy without retardation, 
whereas the full line represents the result with retardation. 

FIG. 4. The free energy vs film thickness for interfacial melting of a methane film on MgO. 
Details as in Fig. 3. 
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TABLES 



TABLE I. Parameters used for an oscillator model of dielectric response for the solid methane 
using data from Ref. [22]. 



ej (eV) fi {eV) gi (eV) 

10.0 10.0 1.8 
12.7 44.0 3.6 
14.4 64.0 3.5 
22.3 10.0 2.4 

23.1 23.0 G.8 



TABLE II. Parameters used for an oscillator model of dielectric response for the liquid methane 
using data from Ref. [22]. 



ej (eV) h {eV^) g^ (eV) 

10.0 8.0 1.5 

12.8 42.0 3.5 

14.8 62.0 4.0 

22.7 8.0 2.4 

23.8 21.0 6.5 



TABLE III. Parameters used for an oscillator model of dielectric response for solid MgO using 
data from Ref. [23]. 



ej (eV) fi {eV^) gj (eV) 

7.6 3.0 0.1 

10.8 150.0 3.5 

13.2 40.0 1.0 

17.0 130.0 4.0 
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Fig.2 
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